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Measurement of v n coefficients and event plane correlations in 2.76 TeV Pb-Pb 

collisions with ATLAS 



Measurements of flow harmonics v\-ve are presented in a broad range of transverse momentum and centrality for 
sJsnn = 2.76 TeV Pb-Pb collisions with the ATLAS experiment. The fourier coefficients of two-particle correlations 
are shown to factorize as products of the single particle v n for n = 2-6. The factorization breaks for n = 1 due to 
effects of global momentum conservation. The dipolar v\ associated with the initial dipole asymmetry is extracted 
from the two-particle correlations via a two component fit that accounts for momentum conservation. The magnitude 
of the extracted dipolar vi is comparable to V3. Measurements of correlations between harmonic planes of different 
orders, which give additional constraints on the initial geometry and expansion mechanism of the medium produced 
in the heavy ion collisions, are also presented. 
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1. Introduction 

In heavy-ion collisions the produced fireball has spatial anisotropies of many orders: elliptic triangular etc. These 
anisotropies are transferred from position space to momentum space due to the collective expansion of the medium 
or path-length dependent suppression of particles. This results in the azimuthal distribution of the final particle yields 
being modulated by the flow harmonics v n about the different harmonic planes 0„ (U : 



The v n can be measured by the event plane (EP) method by correlating the single-particle azimuthal distributions 
with the harmonic planes <D W . They can also be measured by the two-particle correlation (2PC) method where the 
particle-pair distribution in relative azimuthal angle Acp = cp a - cp^ are measured (the subscripts a and b label the two 
particles, commonly called trigger particle and partner particle). The 2PC correlation function can be expanded in a 
Fourier series in A0 as: 
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Figure 1: A sample of V2-ve values measured by the event plane method. The first panel shows the centrality dependence for the v n . The second 
panel shows the pj dependence for (20-30)% centrality and the last panel shows the rj dependence. 



If the correlations are dominated by single particle anisotropics (Eq. [I]), then the Fourier coefficients v„ ?w are equal to 
the product of the individual single particle v n : 

V n ,n(PT> Pj) = V n (j>j)v n (p\) (3) 

Using Eq. [3] one can obtain the v n from the 2PC. The above relation is violated if non-flow effects, such as jets and 
resonance decays are large. Thus the non-flow effects must be suppressed before extracting the flow harmonics from 
the correlations. 

The flow harmonics are important observables as they contain information about the initial geometry and transport 
properties of the medium El [3). A better understanding of the v n can also explain the origin of the ridge, an elongated 
structure along At/ at Acp ~ |4] and the so-called "mach-cone", a double-hump structure on the away-side (5 ] seen 
in the 2PC. These were initially interpreted as the response of the medium to the energy deposited by the quenched 
jets. However, recent studies 1 6 ] have shown that higher-order flow harmonics can contribute to these structures. 

Another set of observables closely related to the v n are the correlations between the event plane angles 0„. These 
correlations can be produced due to correlations between eccentricities of different orders in the initial geometry or 
they can develop during the dynamical expansion of the produced matter. Thus these measurements provide additional 
constraints on initial geometry and the expansion mechanism of the medium (7] [8). 

The v n results presented are for charged hadrons reconstructed in the ATLAS inner detector covering \q\ < 2.5. 
The ®„ planes for the V2-V6 measurements via the event plane method were determined using the forward calorimeter 
covering I77I e (3.3,4.9). For the event plane correlation analysis, the entire EM calorimetry covering 77 e (-4.9,4.9) 
was used. All measurements were done using 8 jib~ l of Minimum Bias Pb-Pb data at of 2.76 TeV. Details of 

the results presented here are published in [10 ] and i fTTll . 

Figure [T] summarizes the results for V2 to V6 measured by the event plane method. Qualitatively the following 
features are seen: the v n coefficients rise and fall with centrality with V2 having considerable centrality dependence 
(being driven by the average geometry). Typically V2 is much larger than the other harmonics, but in most central 
collisions V3, V4 and even V5 can become larger than V2. The v n coefficients rise and fall with pj, this is interpreted 
as the driving mechanism behind the v n changing from collective expansion at low pj to path-length dependent 
suppression at high p^. The v n coefficients are approximately boost invariant showing that they represent non-local 
correlations. 



2. Two-Particle Correlations 

Figure [2] shows the procedure used for obtaining the v n from the 2PC. Panel a) shows the two-dimensional correla- 
tion function in At] - A(p for p\,p\ e (2, 3) GeV. Such correlations, where both trigger and partner particles are in the 
same pj range, are termed as fixed- pj correlations. Long range structures, namely the ridge at Acp = and the double 



/Nuclear Physics A 00 (2012 ) 7-[| 3 




Figure 2: The steps involved in the extraction of the v n for 2-3 GeV fixed-px correlation. Panel a): the two-dimentional correlation function (shown 
for IA77I < 4.75 to reduce fluctuations near the edge). Panel b): The one-dimentional A0 correlation function for 2 < IA77I < 5 (black points) overlaid 
with its individual Fourier components and their sum (thick black line). Panel c): v n vs. n (The bars and bands indicate statistical and systematic 
uncertainties respectively). The lower sub panel in c) shows v n for n > 7 but in a linear scale. 



hump at A0 = n are seen along the Ar] axis. The narrow peak at (A77, A0) ~ (0, 0) is due to jets and other short range 
correlations and is removed by applying a IA77I > 2 cut. Panel b) shows the one-dimensional correlation function in A(p 
which only contains contribution from the long range structures (in Ar]). The ridge and the double hump are clearly 
visible in the 1-D correlation. As the trigger and partner /? T ranges are chosen to be the same, Eq. [3]reduces to: 

v^04,/4) = v„04) 2 (4) 

which is used to obtain the v„(/?t). In Fig. [2j the v n are plotted up to n = 15, but the analyis is limited to n = 6 as for 
higher n the systematic and statistical errors are too large for the v n values to have any significance. 

The 2PC results can be used to check where the factorization (Eq. [3} breaks. If collective flow dominates the 2PC, 
then the near- side peak must be larger than the away-side peak. From Fig. [3] it is seen that this is roughly true up to 
50% centrality (for p\,p\ £ (3,4) GeV) beyond which the away-side peak becomes larger indicating a break-down of 
the v nrn factorization. 

In Fig. [4] the evolution of the correlations is shown for (0 - 10)% central collisions. At low pj the correlation 
is driven by flow and the near- side peak is larger than the away-side. However, for p\,p\ > 6 GeV, the 2PCs are 
dominated by the away-side jet peak. Thus it is clear that at high p^ (> 6 GeV) and in peripheral collisions (> 50% 
centrality), non-flow effects become important and Eq. [3]is not expected to hold. 



3. Comparison of v n obtained from 2PC and EP methods 

The left panels of Fig. [5] compare the V2 from fixed-/?T 2PC method with those from the EP method for (0 - 10)% 
central collisions. The two methods agree within 5% to 15% for V2 for pj < 4 GeV. Deviations are observed for 
p T > 4 GeV, presumably due to contributions from the away-side jet. 

The v n values obtained from fixed-pr correlations can be cross-checked using correlations where the trigger and 
partner have different p\ and p^ values. Such correlations are termed as mixed- pj correlations. The v n (pj) of the 
partner can be obtained as: 

Equation [5] can be checked by measuring the same v n (p\) for different v n (pj). This is illustrated for V2 in the right 
panel of Fig. p] It is seen that factorization of v n ^ n works well for n = 2 and the values of V2(p\) are reasonably 
independent of p\. Further it is seen that when p* is below 3 GeV, the agreement of the 2PC values with the EP 
method extend out to much higher p^. This shows that as long as the reference p\ is low, the factorization relation 
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Figure 3: Centrality dependence of A0 correlations for p\,p\ e (3,4)GeV and A77 e (2,5). The superimposed solid lines (thick-dashed lines) 
indicate contributions from individual vn,n components (sum of the first six components). 




Figure 4: Fixed-px correlation functions in the (0-10)% centrality interval for several pj ranges with 2 < \At]\ < 5. 




p; [GeV] p> [GeV] 



Figure 5: Left panels: comparision between the fixed-px ^i(pt) (or V3(px)) and EP V2(px) (or V3(px)- Right panels: Comparison of V2(p\) (y^pr)) 
obtained for four different reference p^ ranges (0.5-1, 1-2, 2-3, 3-4 GeV) with the EP values. The error bars indicate the statistical uncertainties 
only. The dashed lines in the ratio plots indicate a ±10% band to guide eye. All results are for the (0-10)% centrality interval. 
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ror bars indicate the statistical un- 
certainties only. 



is valid and hence the v n can be measured to high pj via the 2PC method. This factorization between hard and soft 
particles is expected since the high pj particles, due to path-length dependent energy loss, are correlated with the 
same geometry that drives the collective expansion at low pj 1 12]. 

In Fig. [6] the validity of the factorization is checked for all harmonics as a function of IA77I for four different p T , p T 
combinations. For harmonics V2 to v$ and \Ar]\ > 1.0, all trigger-partner combinations give nearly identical values of 
v n (p T ) even though the trigger v n values vary over a large range. This validates the factorization for n = 2 - 6. 



4. Extracting dipolar v\ from v\ 9 \ 

The first panel of Fig. [6] shows that for vij the factorization does not hold for any value of the Arj gap. This is 
because the vij is strongly influenced by global momentum conservation (GMC), which modifies the factorization 
forv u to O: 

vuO^t^V) = vi(p*,i7>i(/4,i7 a ) - "S^r (6) 

M</4> 

where, M is the multiplicity in the event. The GMC term in the above equation is the leading order approximation 
for momentum conservation and is important at high p T and when the multiplicity is low. The v\(p^,7f) can be 
decomposed into rapidity-even and rapidity-odd components. The rapidity-odd component of vi is due to sideward 
deflection of the colliding nuclei and is small at mid-rapidity (less than 0.005 for I77I < 2), hence its contribution to vi ? i 
is small ( < 2.5 x 10" 5 ). The rapidity-even component comes from a dipole asymmetry due to fluctuations in in the 
initial geometry and like the other v n is expected to be large and have a weak 77 dependence 1 14 ]. In this case, Eq. [6] 
simplifies to: 

vu(Pt'/>t) = vi04)vi04) - cp*p\ (7) 
where c = M ^ 2 - . The v\(pj) are obtained as fit parameters by fitting vij(^, p\) according to the above equation with 

c as an additional fit parameter. Figure [7] shows the results of the fit for (0-5)% central collisions. The lower panels in 
Fig. [7] show the difference between the fit and the data indicating that the two component fit works very well. 

Figure [8] shows the dipolar vi(/?t) values obtained from the two component fit with the V2 and V3 values. A large 
dipolar v\ is observed and is comparable to V3, indicating a significant dipolar asymmetry in the initial geometry. The 
Pj dependence of the dipolar vi is similar to other harmonics: it increases up to intermediate pj then decreases, which 
is expected as its origin is the same as the other harmonics. The negative dipolar v\ seen at low p? is expected from 
hydro calculations 031 • 




Figure 7: Global fit to the vij data for the (0-5)% centrality interval. The fit is performed simultaneously over all vi,i data points in a given 
centrality interval, which are organized as a function of for various ranges (indicated at the top of each panel), with shaded bars indicating 
the correlated systematic uncertainties. The fit function and the associated systematic uncertainties are indicated by the thick-solid lines and 
surrounding dashed lines, respectively. The dot-dashed lines intercepting each dataset (with negative slope) indicate estimated contributions from 
the momentum conservation component. The lower part of each panel shows the difference between the data and fit (solid points), as well as the 
systematic uncertainties of the fit (dashed lines) 



The agreement between the 2PC and EP methods implies that the structures of the two-particle correlation at low 
Pt and large At] mainly reflect collective flow. This is verified explicitly by reconstructing the correlation function as: 

C reco (A0) = N^ c ll + 2vlfcos(A(f>) + 2 ^ vf vf cos(nA^)\ (8) 

\ n=2 > 

where Nq FC and are the average of the correlation function and first harmonic coefficient from the 2PC analysis, 
and the remaining coefficients are calculated from the v n measured from the event plane method. Figure [9] compares 
two measured 2PCs to the corresponding reconstructed correlations (Eq. [8]) showing excellent agreement between the 
two. The 2PC in the right panel of Fig. [9] has a large vi,i component which contributes to the double hump, however 
Fig. [8] shows that a large fraction of the Vi,i at this pj comes from the dipolar Vi. This demonstrates that the ridge and 
cone are a manifestation of single particle v\-v^. 



5. Event plane correlations 

Further insight into initial geometry can be obtained by studying correlations between the O w . The correlations 
between two angles O w and O m are described by the differential distribution ) , where k is the lowest common 

multiple of n and m. This distribution can be expanded as a Fourier series as GO [8): 

Jl^l \ = 1 + 2I >1 y w <*>*UW>n " ®m)) (9) 

dk(<3> n - O m ) 3 

where the Fourier coefficients V^ m quantify the strength of the correlations. The measured Fourier coefficients need 
to be corrected to account for the detector resolution effects ifTTTl . 

The two-plane correlations are summarized in Fig. [TOj The first three panels show the j = 1,2 and 3 moments of 
the 4(®2 - ^4) correlation as a function of the number of participating nuclueons (A^ par t). The correlations are small in 
most central collisions and increase almost linearly with decreasing (N part ) becoming fairly large in mid-central and 
peripheral collisions showing a strong correlation between O2 and O4. 
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Figure 9: Measured correlation functions compared 
with those reconstructed from V2-ve measured with 
the EP method and vij with the 2PC method in (0- 
1)% centrality interval for two pj ranges. The con- 
tributions from n=l, n=3,5 and n=2,4,6 are shown 
separately. The error bars indicate the statistical un- 
certainties. 



The fourth panel shows the j = 1 moment of the correlation between O2 and ® 3 . The correlations are very weak, 
consistent with zero in the most central collisions and increase to roughly 2% in peripheral collisions. The next two 
panels show the j = 1 moments of the O2 - ®6 an d ^3 - ®6 correlations. The O2 and O3 are weakly correlated with 
each other, but they individually are strongly correlated with $5. Also the two correlations show completely different 
centrality dependance : O2 - ®6 correlation increases almont linearly with decreasing (N part ) (i-e. from central to 
peripheral collisions) while O3 - ®6 gradually decreases. The last two panels of Fig. 10 show the O3 - O4 and O2 - O5 
correlations respectively, which are found to be weak (less than a few %) and consistent with zero. It is possible to 
generalize the two-plane correlations to correlations between three or more planes. A detailed analysis of these event 
plane correlations is presented in a separate talk in this conference by J. Jia fT5l . 



6. Summary 

Measurements of the flow harmonics v\-ve over a large pj, 77 and centrality range are presented. The 2PC v n , n are 
shown to factorize into products of single particle v n for n = 2 - 6 in central and mid-central collisions as long as one 
particle has low p? (< 3 GeV). The factorization breaks for vij due to the contribution from momentum conservation 
effects. Dipolar v\ is extracted from the vi,i data via a two component fit. The extracted dipolar vi is comparable 
to V3, indicating significant dipole asymmetry in the initial geometry. It is shown that the features in two-particle 
correlations for |At/| > 2 at low and intermediate p? (< 3 GeV) such as the ridge and double-hump are accounted for 
by the vi-v&. 

The v n can be thought of diagonal components of a Flow Matrix. Studying the two and three plane correlations 
gives access to the off diagonal entries and beyond. A detailed set of event plane correlation measurements are 
presented. 
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Figure 10: Two-plane correlators (cos jl(O n -O m » as a function of <A^ pa rt>- The middle two panels in the top row have j = 2 and j = 3, respectively, 
while all other panels have 7=1. The error bars and shaded bands indicate the statistical and systematic uncertainties respectively. 
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